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We theoretically study homogeneously trapped atomic Bose-Einstein condensates where all three 
momentnm components conple to a pseudo-spin-1/2 degree of freedom. Tuning the anisotropies of 
spin-orbit coupling and the spin-dependent interactions is shown to provide access to a rich phase 
diagram with a tetracritical point, first-order phase transitions, and multiple metastable phases of 
stripe and plane-wave character. The elementary excitation spectrum of the axial plane-wave phase 
features an anisotropic roton feature and can be used to probe the phase diagram. In addition 
to providing a versatile laboratory for studying fundamental concepts in statistical physics, the 
emergence of metastable phases creates new opportunities for observing false-vacuum decay and 
bubble nucleation in ultra-cold-atom experiments. 
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The possibility to create artificial gauge fields in neu¬ 
tral ultra-cold atom systems 0, i has drastically ex¬ 
panded the array of possibilities for highly controlled 
experimental simulation of quantum many-particle sys¬ 
tems ii- In particular, it has become possible to ex¬ 
plore effects associated with spin-orbit coupling (SOC) 0 
that give rise to intriguiM phenomena such as the 
quantum spin Hall effect (5l-l8|, new materials classes 
such as topological insulators and superconductors 0- 
and exotic qu asip article excitations such as Majo- 
rana fermions |lll - ll3l| . In bosonic-atom systems, the 
presence of SOC was found to generate novel ground 
states that have no known analogues in conventional 
solid-state materials 14-1^. Intense theoretical atten¬ 
tion has focused on the many-body physics of spin-orbit- 


coupled Bose-atom systems in free space |17H23l| and in 
harmonic traps [2J-l30|. To date, only a special type 
of SOC involving a single Cartesian component of the 
atoms’ momentum has been realized in the lab 3ll-34|. 
However, several proposals exist for creating a Rashba- 
type SOC in higher dimensions 35l-l37j|. Experimen¬ 
tal progress is spurred by predictions of an exotic half¬ 
vortex phase [2^, [26j and a striped superfluid phase 15 1 
is systems with two-dimensional Rashba SOC. The three- 
dimensional analogue of Rashba SOC is also interesting 
because it is expected to stabilize a long-sought Skyrmion 
mode in the gro und state of trapped Bose-Einstein con¬ 
densates 271, l30| . The intriguing possibility to simulate 
the so-called Weyl SOC has also been suggested [35 1. 


Despite the rapid pace of theoretical and experimen¬ 
tal studies in the field of spin-orbit-coupled atom gases, 
the physical properties of an extended interacting Bose 


system in the presence of three-dimensional (3D) SOC 
have not yet been considered in any detail. This clearly 
presents a crucial gap in our basic understanding, as 
the extended system’s behavior constitutes an important 
benchmark for identifying effects associated solely with 
trapping potentials. Furthermore, real experimental sys¬ 
tems can be designed with a flat-bottom potential to ap¬ 
proximate the extended system and give access to the in¬ 
triguing physics demonstrated by our present study [s^ . 
In particular, we show that 3D SOC in an interacting 
Bose gas leads to a highly nontrivial phase diagram fea¬ 
turing a tetracritical point, first-order phase transitions, 
and emergent metastable phases - none of which have 
been seen in systems with lower-dimensional SOC and/or 
in the presence of a trapping potential. Thus this system 
provides opportunities to study ramifications of multi¬ 
criticality |39| and metastabili^ including false-vacuum 
decay and bubble nucleation |4Cll443l| . in ultra-cold-atom 
experiments. We have also studied the spectrum of ele¬ 
mentary excitations and find it to be useful for probing 
the multitude of phases and phase transitions. 

The model. — We consider a 3D homogeneous inter¬ 
acting two-component Bose gas subject to cylindrically 
symmetric spin-orbit coupling, described by the Hamil¬ 
tonian H = Hq + Hi, with 


Ho= d^r 
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Here ^'(r) = is a two-component spinor field, 

Ua = ipti’cr is the density for component cr G is 

the atomic mass, aj (with j = x, y, z) denote the Pauli 
matrices and pj = —ihVj are the Cartesian components 
of the single-atom momentum operator p. 

The parameter A measures the SOC strength involv¬ 
ing the momentum pj^ in the xy plane, and the dimen¬ 
sionless number 7 describes the anisotropy of SOC for 
the momentum component parallel to the z direction. 
Note that the limit 7 = 0 is the unitary equivalent of 
the conventional Rashba SOC [l|, 7 = 1 realises the so- 
called Weyl SOC [s^, and a situation corresponding to 
the experimentally created Rashba-type SOC |8ll434l | is 
obtained when 7 —>■ 00 (with A 7 finite). Fundamentally, 
the parameter 7 could be tuned by a sequence of pulsed 
inhomogeneous magnetic fields [^. Also, the strength g 
(Sti) interactions between same-spin (opposite-spin) 
components can be varied using an appropriate Fesh- 
bach resonance In the special case when 7 = 1 

and g = g-i-j,, the Hamiltonian H is symmetric with re¬ 
spect to simultaneous rotations of the internal pseudo¬ 
spin- 1/2 degree of freedom and the particle momentum. 
Throughout the rest of the paper, we use units such that 
h = Ub = 2 to = 1 . 

Diagonalization of Hq yields the two-branch single¬ 
particle energy spectrum E±{p) = (p± Arp/2)^ — A^Tp/4 
as a function of 3D momentum p. Using spherical coor¬ 
dinates, p = (psin 0 p cos(pp,psin 0 p sin(/ 9 p,pcos 0 p), we 

have Tp = 7 ^ cos^ 6p + sin^ 9p, and the eigenspinors 

are given by 


± sin 0D 


.p±(p) = 




\/A> 7 
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\/ 2 Tp 
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The lowest-energy state for a given propagation direction 
parameterized by dq and ipq is from the ” branch and 
occurs at the momentum q satisfying q = ATq/ 2 . 

The phase diagram. — To determine the ground state 
of the interac ting system, as it is routinely done in the 
literature [II, [ 2 ^, we assume that the system has con¬ 
densed into a coherent superposition of two plane-wave 
states with momenta ±q having magnitude q = ATq/ 2 . 
Thus the condensate wave function has the form $0 = 
C+ 'I'_(q) -I- C- '!'_(—q), with coefficients C± that will 
be determined by a variational procedure. The condi¬ 
tion no = lU+p -I- lU-p, with uq being the particle 
number density, suggests the parameterization \C-\'^ = 
noCOS^(a/2) and jC-p = n-osin^ (a/2), with a G [ 0 , 7 r]. 
Introducing the dimensionless nonlinear-coupling param¬ 
eter g = ((/-[-j, — g)no/X^, we find the variational ground- 
state energy density Eg = ($o|i?|<i>o)/V — griQ given by 


A^no 
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FIG. 1. (color online) The ground-state energy density Eg 
(measured in terms of A^no) plotted as a function of the vari¬ 
ational parameters dq and a, for particular values of quan¬ 
tities characterising anisotropy of nonlinear interactions (g) 
and spin-orbit coupling ( 7 ): (a) g = 1.0 and 7 ^ = 0 . 8 ; (b) 
g = 0.25 and 7 ^ = 0.5; (c) g = —0.1 and 7 ^ = 1.8; and (d) 
g = —0.25 and 7 ^ = 0.5. The global minimum in each panel 
corresponds to a true ground state, while the existence of lo¬ 
cal minima on the energy landscape signifies the emergence 
of metastable phases. 


We first consider the familiar case of two-dimensional 
SOC by setting 7 = 0 . Minimization of Eg with respect 
to 9q and a then yields sin^ 0 q = 1 and sin^ a = 1 ( 0 ) 
for g > 0 (g < 0). The first condition implies that the 
condensation momentum is pinned in the xy plane, and 
the latter condition yields the stripe phase for g > 0 
(1(7 +1 = |(7_|, i.e. condensation in a coherent superposi¬ 
tion of the opposite-momentum states) or the plane wave 
phase for g < 0 (either |(7+| = 0 or \C-\ = 0, i.e. con¬ 
densation at only one momentum eigenstate). We thus 
reproduce the ground-state structure of the conventional 
two-dimensional-Rashba SOC case 0 . 

Setting 7 yf: 0 unpins the condensation momentum 
from the xy plane, making it possible to condense into 
a state whose momentum has a finite z component. In 
the following, we will term such condensation as “Polar”, 
while condensation into a momentum that lies in the xy 
plane will be called “Axial”. As each of these cases can 
support a stripe (SP) or plane-wave (PW) condensate, 
depending on the interaction strength g, we have four dis¬ 
tinct possible phases: PW-Axial, SP-Axial, PW-Polar, 
and SP-Polar. Examination of the variational ground- 
state energy landscape shows that each of the four phases 
is found to be either a true ground state or a metastable 
state, depending on the values of g and 7 . See Fig. [T] for 
pertinent examples. The stability of the phases can be 
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FIG. 2. (color online) Phase diagram controlled by vary¬ 
ing the parameter g that measures anisotropy of spin- 
dependent interaction strengths and the quantity 7 ^ related to 
anisotropy of spin-orbit coupling. Four possible phases exist 
- PW-Polar, SP-Polar, PW-Axial, and SP-Axial - that inter¬ 
sect at the tetracritical point Qc- For Polar (Axial) phases, 
the condensate momentum has a finite (vanishing) z compo¬ 
nent. In the PW (SP) phases, condensation occurs into a 
single (a superposition of two) plane-wave state(s). 
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FIG. 3. (color online) The phase diagram of metastable 
states controlled by g and 7 ^. White regions indicate pa¬ 
rameter combinations for which there no dynamically sta¬ 
ble metastable phases exist, i.e., where the Hessian matrix 
in Eq. O is not positive-definite. 


ensured by the positivity of the Hessian matrix 


He = 


a^E„ dE„ aE„ 

as^ aa 

dE„ aE„ a^E„ 


(4) 


dOa da 


da'^ 


The true ground state phase diagram spanned by tun¬ 
ing parameters g and 7 ^ is shown in Fig. [S] A tetra¬ 
critical point Qc connecting the four possible phases 
emerges when g = 0 and 7 = 1 . At this high-symmetry 
point, the system is invariant with respect to simultane¬ 
ous SU(2) spin rotation and rotation of the momentum 
of the atoms. The observed behavior at Qc in our sys¬ 
tem contrasts with that exhibited in the presence of a 
tight harmonic trapping potential where Skyrmion tex¬ 
tures are stabilized in the ground state [27j |. 

The lowest dynamically stable metastable states are 
shown in Fig. [31 They gradually disappear as the pa¬ 
rameter 7 ^ approaches 0 , in the sense that local minima 
in Fig. |T] cease to exist at this point. This means that 
metastable phases literally emerge in Bose-Einstein con¬ 
densates with 3D SOC only. The presence of metastable 
phases along with the true ground states creates the op¬ 
portunity to simulate false-vacuum decay. Proposed by 
Sydney Coleman for modeling phase transitions in the 
universe [4l|, decay from a false vacuum into a true one 
plays a key role in numerous physical contexts. For ex¬ 
ample, it occurs in a superheated liquid, where the false 
vacuum is the liquid state, while the true one is gaseous 
[^ . Thermodynamic fluctuations trigger the continu¬ 
ous appearance of vapor bubbles in the liquid. Eventu¬ 
ally growing bubbles swallow the entire system. More 
speculative manifestations of the phenomena exist also 
in modern cosmology 4^, 461. Due to the its high tun- 
ability, our system provides an easy route toward testing 


the false-vacuum quantum decay. The system can be pre¬ 
pared initially in one of the metastable phases of Fig. [S] 
Quantum fluctuations are then expected to trigger quan¬ 
tum decay accompanied with nucleation of bubbles of one 
of the lower-lying true ground states. 

Elementary excitations. — The phases and phase tran¬ 
sitions in our system can be probed by studying the spec¬ 
trum of elementary excitations, e.g., by using Bragg spec¬ 
troscopy Here we consider the elementary exci¬ 

tations around the PW-Axial ground state; partly mo¬ 
tivated by the fact that, for the case of one-dimensional 
SOC, interesting roton-like modes were found Phys¬ 
ically, the roton mode signals a system’s tendency to un¬ 
dergo a first-order phase transition to a supersolid when 
the roton gap closes [13, [Hj , and it is usually the con¬ 
sequence of strong correlations in the system due to the 
interplay of SOC and interactions. Our aim is to show 
that these features persist also in the case of 3D SOC and 
that it probes the rich phase diagram obtained above. 

The PW-Axial phase has one condensation momen¬ 
tum lying in the xy plane. Without loss of gener¬ 
ality, we choose the condensate momentum to be = 
■|(—1,0,0). Within the framework of imaginary-time 
functional integration, the partition function of the sys¬ 
tem reads Z, = j T] exp (—'F]) with the 

action Tf] = dr)/ (Pr Y,a + H - fiN], 

where /3 = 1/T is the inverse temperature, and fi is 
the chemical potential introduced to fix the total par¬ 
ticle number. The Bose field is split into the mean-field 
and fluctuating parts, = $0CT|q^K + 
expand the action of the system up to the quadratic order 
in fluctuating fields obtaining an effective action SeS — 

So + Sg. Here S'o = V (“X “ ^)^ 0 rT + (fl + gn)'<^la 
is the mean-field contribution, while Sg = 74>lC/“^<Fq 
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is the fluctuating contribution with a vector field 
$q = IS the inverse 

Green’s function of the elementary excitations defined as 

f-iwn + e+ i?q gno g\ino \ 

g-i ^ Rq -iwn + e” g^ino gno 

~ gno guno iwn + elq Rlq 

V gnno gno R^q iWn + eZqJ 

(5) 

where Cq = + ^ + A(±7g2 — Qx) + gno and i?q = 

gnno + A(-f + qx- iqy)- 

The spectrum of the elementary excitations is deter¬ 
mined from the poles of the Green’s function. There are 
two branches of excitations found as illustrated in Fig. 2) 
We choose parameters to probe the PW-Axial ground 
state and show the results in the first row in Fig. 21 The 
lower branch of the spectrum exhibits a typical linear 
Bogoliubov slope at low momenta, followed by a roton 
and maxon features at higher momenta. This structure 
of the spectrum persists whenever the momentum q has 
components along the axis of the condensation, in this 
case the cc-axis. The roton-maxon feature is absent along 
any direction that is perpendicular to the direction of 
condensation momentum. There is a conical intersection 
at around = 1, reflecting the time-reversal symmetry 
of the system. It can be lifted by a Zeeman-like field in 
the Hamiltonian, in which case the lower branch will be¬ 
come separated from the upper one. Our purpose here 
is to study how the spectrum changes when we drive the 
system across the phase diagram. 



FIG. 4. (color online) Dispersion of low-lying elementary ex¬ 
citations u)± along Qx (left panels), Qy (centre panels), and qz 
(right panels). First row: g = —0.5 and 7 ^ = 0.5 (SP-Polar 
phase is stable), a phonon-maxon-roton feature is seen along 
the X direction. Second row: g = 0 and 7 ^ = 0.5 (boundary 
between PW-Axial and SP-Axial phases), the roton minimum 
goes soft. Third row: g = —0.2 and 7 ^ = 1.5 (metastability 
of the PW-Axial phase is broken), imaginary parts appear in 
qz dispersion. We have set gno/X^ = 0.25 in all panels. 


When the system is driven close to the boundary with 
the SP-Axial phase, the roton minimum becomes soft as 
shown in the second row of Fig. 2) This signals the in¬ 
stability of the system toward the striped order, whose 
spatial modulation is set by the momentum at which 
the gap closes. There are no metastable phases at this 
boundary, as shown in Fig. [3] This phase transition is 
of second order. On the other hand, nothing dramatic 
occurs when we move close to the boundary with the 
SP-Polar phase. The roton gap does not close, and the 
spectrum of the PW-Axial phase does not show any spe¬ 
cific feature at the phase boundary. This is due to the 
presence of metastable phases. The PW-axial phase be¬ 
comes metastable when we cross the line 7 ^ = 1 at fixed 
g < 0 from below, as shown in Fig. [3] Therefore, this 
phase transition is of the first order. We need to drive 
the system much further to see changes in the excita¬ 
tion spectrum, namely until the point when metastabil¬ 
ity breaks down as it is shown in the third row of Fig. 21 
The spectrum of the SP phase is qualitatively different 
from the PW-axial phase and features a double-gapless 
band structure due to spontaneously broken translations 
symmetry as shown in Ref. Therefore, measuring 

the excitation spectrum can be used as a probe of the 
rich phase diagram in the presence of 3D SOC. 


It is interesting to note that at the tetracritical point 
Qc, we find also two gapless Goldstone modes, result¬ 
ing however from spontaneous breaking of spin-rotation 
symmetry and U(l) gauge symmetry. Such modes are ex¬ 
pected to remove the four-fold degeneracy at Qc found at 
the mean-held level, leading to a unique ground state via 
the so called “order from disorder” mechanism 26, 531. 


Experimental relevance. — For a trapped Bose gas in 
the presence of Weyl SOG and weak inter-particle inter¬ 
action, one expects that the ground state is a Skyrmion, 
which is a superposition of few lowest Landau levels . 
Our pre dictions should apply for hat bottom traps as in 
Ref. [38l |. In addition, we may expect the main features 
of the presented phase diagram be present in harmonic 
traps with sufficiently strong nonlinearity [s^. 
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